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ABSTRACT 

Aims. We have obtained high-resolution spectroscopy, optical interferometry, and long-term broad band photometry of the ellipsoidal 
primary of the RS CVn-type binary system f And. These observations are used to obtain fundamental stellar parameters and to study 
surface structures and their temporal evolution. 

Methods. Surface temperature maps of the stellar surface were obtained from high-resolution spectra with Doppler imaging tech- 
niques. These spectra were also used to investigate the chromospheric activity using the Ha line and to correlate it with the pho- 
tospheric activity. The possible cyclicity in the spot activity was investigated from the long-term broad band photometry. Optical 
interferometry was obtained during the same time period as the high-resolution spectra. These observations were used to derive the 
size and fundamental parameters of £ And. 

Results. Based on the optical interferometry the apparent limb darkened diameter of f And is 2.55 ± 0.09 mas using a uniform disk 
fit. The expected ~4% maximum difference between the long and short axes of the ellipsoidal stellar surface cannot be confirmed 
from the current data which have 4% errors. The Hipparcos distance and the limb-darkened diameter obtained with a uniform disk fit 
give stellar radius of 15.9 + 0.8R o , and combined with bolometric luminosity, it implies an effective temperature of 4665 ± 140 K. The 
temperature maps obtained from Doppler imaging show a strong belt of equatorial spots and hints of a cool polar cap. The equatorial 
spots show a concentration around the phase 0.75, i.e., 0.25 in phase from the secondary, and another concentration spans the phases 
0.0-0.4. This spot configuration is reminiscent of the one seen in the earlier published temperature maps of f And. Investigation of 
the Ha line reveals both prominences and cool clouds in the chromosphere. These features do not seem to have a clearly preferred 
location in the binary reference frame, nor are they strongly associated with the cool photospheric spots. The investigation of the 
long-term photometry spanning 12 years shows hints of a spot activity cycle, which is also implied by the Doppler images, but the 
cycle length cannot be reliably determined from the current data. 

Key words, stars: activity - stars: chromospheres - stars: fundamental parameters - stars: individual: f And stars: starspots 
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' 1. Introduction magnetic activity can only be produced when differential rota- 

f— 3 ' „ , , , , . , ,. , tion is present (Riidiger et al. l2003l . 
^— "j Because or the enhanced dynamo action in stars with thick, tur- 

. . bulent outer-convection zones, rapidly rotating cool stars, both In this work f Andromedae, a long-period (17.8 day), 

^ ■ evolved and young, exhibit significantly stronger magnetic ac- single-lined spectroscopic RS CVn binary (Campbell 1191 It 

tivity than is seen in the Sun. This activity means that the spots Cannon 1191 5b , is investigated in detail. In this system the pri- 

are also much larger than the spots observed in the Sun. The mary is of spectral type Kl III, and the unseen companion pos- 

■ largest starspot recovered with Doppler imaging is on the active sibly of type F (Strassmeier et al. 119931) . The primary fills ap- 

RS CVn-type binary HD 12545 which, in January 1998, had a proximately 80% of its Roche lobe, so it has a non-spherical 

spot that extended approximately 12x20 solar radii (Strassmeier shape. The estimated ellipticity gives a ~4% difference between 

1999 ). The lifetime of the large starspots/spot groups can also the long and short axes of the ellipsoid (Kovari et al. 120071 from 

be much longer than that of the sunspots, even years instead here on Paper 1). The mean angular diameter of £ And has been 

of weeks for sunspots (e.g., Rice & Strassmeier 1 19961 Hussain derived to 2.72 + 0.036 mas using spectro-photometry (Cohen et 

2002 ). The most typical dynamo signature is the presence of an al. 1 19991 ) 
activity cycle. Cyclic changes in the level of magnetic activity 

are well documented for the Sun, as well as for many solar-type An earller detalled Do PP ler lma g ln 8 stud y < Pa P er ! > revealed 
stars (see, e.g., Olah et al. [20091). It is also interesting that, ac- that the s P ots on the surface of f And have a temperature con- 
cording to theoretical calculations, cyclic variations in the stellar ttast of approximately 1000 K and that they occur on a wide 

latitude range from the equator to an asymmetric polar cap. The 

* Based on observations collected at the European Southern strength of the features changed with time, with the polar cap 

Observatory, Chile (Prg. 081.D-0216(AB)); and with the Vienna auto- dominating the beginning of the two-month observing period in 

made photometric telescopes Wolfgang and Amadeus, Arizona, USA 1996/97, while the activity during the second half was domi- 
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nated by medium-to-high latitude features. Also, the investiga- 
tion revealed a weak solar-type differential rotation. 

Here, results from Doppler imaging, optical interferometry, 
and long-term photometry of £ And are presented. We discuss 
the reduction of the interferometric data and the obtained funda- 
mental stellar parameters. The high-resolution spectra are used 
with Doppler-imaging techniques to obtain a surface temper- 
ature map. This surface map is compared to the earlier pub- 
lished temperature maps and also with the chromospheric activ- 
ity based on observations of the Ha line. Finally, the long-term 
broad band photometry is used to study the temporal evolution 
of the spottedness, hence the possible spot cycles. 

2. Observations 

Simultaneous observations were carried out at the European 
Southern Observatory with UVES (UV- Visual Echelle 
Spectrograph; Dekker et al. 2000) mounted on the 8-m 
Kueyen telescope of the VLT, and the AMBER (Astronomical 
Multi BEam combineR; Petrov et al. 120071) instrument of the 
VLT Interferometer (VLTI). Additionally broad and intermedi- 
ate band photometry in V, Ic and y bands were obtained with 
the automatic photoelectric telescopes Wolfgang and Amadeus 
in Arizona, USA (Strassmeier et al. 119971 Granzer et al. I20011 >. 
For all the photometric observations HD 5516 was used as the 
comparison star. 

All the observations were phased using the same ephemeris 
as in Paper 1 , 

HJD = 2 449 997.223 + 0.017 + 17.769426 + 0.000040 x E, 
referring to the time of the conjunction. 

2. 1 . Optical interferometry 

The AMBER observations were obtained during the second part 
of the nights starting on September 14, 16, and 18, 2008, corre- 
sponding to orbital phases <p - 0.05 (secondary in front), cp = 
0.15 (intermediate case), and <p = 0.27 (secondary to the side), 
respectively. The details of the observations are listed in TableQ] 
only available online. During the night starting September 18 the 
coherence time was very short, so the data quality is lower than 
during the other half nights. For all the observations, AMBER 
was used in the low-resolution mode at /, H, and # passbands, 
giving a resolving power (A/AA) of ~ 35 and recording data be- 
tween about 1. 1-2.5 /im. Only the H and # band data (~ 1.5- 
2.5 /mi) were used for the data analysis. The J band data were of 
poor quality owing to vanishing detected flux. The fringe tracker 
FINITO (Le Bouquin et al. 2008 ) was used for most observa- 
tions. During the night starting September 14, data were also 
taken without the use of FINITO in order to confirm the cali- 
bration of the visibility. The Auxiliary Telescopes (ATs) were 
placed at the stations AO, K0, and Gl, giving ground-baseline 
lengths of 128m (A0-K0) and 90m (A0-G1 and K0-G1). The 
A0-G1 and K0-G1 baselines have the same ground length, but 
differ in position angle by 90 deg. 

In addition to £ And, a circular check star was observed 
every night. For this // Peg was chosen because it is at a sim- 
ilar position on the sky as £ And and it is expected to have 
a similar angular diameter (@ld=2.50 + 0.04 mas; Nordgren et 
al. 1200 II Mozurkewich et al. 12003b . Observations of £ And and 
fi Peg were interleaved with observations of the interferometric 
calibration stars 6»Psc (Kl III, #=1.86, LD =2.OO ± 0.02 mas), 
41Psc (K3III, #=2.43, © LD =1-81 ± 0.02mas), HD7087 



(G9III, #=2.48, LD =l-59 + 0.02 mas), and HD 15694 (K3 III, 
#=2.48, ©ld=1-77 + 0.02 mas). The angular diameters for#Psc 
and 41 Psc are from Borde et al. (2002) and those for HD7087 
and HD 15694 are from Merand et al. ( 120061) . 

2.2. Spectroscopy 

The UVES observations of £ And were carried out during 10 
nights between September 13, 2008 and October 1, 2009. The 
red arm in the standard wavelength setting of 600 nm was used 
with the imageslicer #3. This instrument setup gives a spec- 
tral resolution (A/AA) of 1 10000 and a wavelength coverage of 
5000-7000 A. Each observation consists of three exposures of 
8 seconds that were later combined to one very high signal-to- 
noise ratio (S/N) spectrum. The S/N of combined observations 
was between 586 and 914 around 6400 A. The data were reduced 
using the UVES pipeline. A summary of the spectroscopic ob- 
servations is given in the on-line Table [2] 

3. Reduction and analysis of the interferometric 
data 

3.1. Data reduction 

Raw visibility and closure phase values were computed using 
the latest version of the amdlib data reduction package (version 
2.2) and the yorick interface, both provided by the Jean-Marie 
Mariotti Center (JMMC). The data reduction principles are de- 
scribed in Tatulli et al. (l2007l l. 

Absolute wavelength calibration was performed by correlat- 
ing the raw spectra with a model of the atmospheric transmis- 
sion, resulting in a correction of AA/A = -0.043 in the K-band 
with respect to the original wavelength table (cf. Wittkowski et 
al. I20081 ). For each observation mentioned in TableQ] only some 
of the individual frames were selected for further analysis. Only 
those frames were used that had a flux ratio under 3 between the 
telescopes of the concerned baseline and that had an estimated 
absolute piston of less than 4 /im. Finally, out of these only the 
30% of the frames with the highest fringe signal-to-noise ratio 
were kept. The selected frames were averaged. 

The resulting differential phase and visibility values were 
significantly affected by chromatic piston effects caused by the 
dispersion of the air (cf. Millour et al. 2008; Le Bouquin et al. 
2009). This effect was relatively strong for our data because of 
the combination of long baselines and large airmasses. We used 
the measured differential phase to estimate the amount of chro- 
matic piston 6 using 5=1 /27r d<pldcr, where <p is the differential 
phase and <x the wavenumber. The loss of the squared visibil- 
ity amplitude p was estimated using formula (1) of Millour et 
al. (f2008 ). The averaged visibility data were compensated using 
the estimated p. Millour et al. note that 5 is the absolute piston 
value relative to the white light fringe. The absolute piston also 
includes the frame-by-frame piston that is estimated by the regu- 
lar AMBER data reduction. This quantity is determined with re- 
spect to the pixel-to-visibility matrix (P2VM) reference, which 
can have an offset to the white light fringe. We selected frames 
with an estimated piston of less than 4 /mi, and verified that the 
piston of our P2 VM measurements is less than 2 /mi in the H 
band and less than 4 //m in the # band. In total, we assumed an 
error of the piston estimate of 5 /mi and propagated it to the final 
visibility amplitude. We also used an alternative compensation 
of the loss of the squared visibility amplitude p that was based on 
a parametrization of the calibrator star data as a function of opti- 
cal path difference, i.e., an estimate that does not depend on the 
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VLTI/AMBER visibility data of f And (top) and of the 
star fi Peg (bottom), compared to uniform disk models. 



Table 3. Uniform disk fit results for the VLTI/AMBER data. 
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± 0.21 mas 
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± 0.06 mas 
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ON 
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±0.11 mas 
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± 0.19 mas 
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ON 
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± 0.12 mas 


2.62 


± 0.16 mas 


#UD 


all data 
all data 


2.48 
2.55 


± 0.09 mas 
± 0.09 mas 


2.43 
2.49 


±0.11 mas 
± 0. 1 1 mas 



cause the data quality is significantly lower on the night starting 
on September 18 than during the two other observing nights. 

During the night of September 14, we obtained data with 
and without the use of the fringe tracker FINITO. The results 
for these two data sets agree well within the errors for both tar- 
gets, and we do not see any systematic calibration effects that 
are caused by the use of FINITO. Deviations between observed 
visibility values and the UD model are mostly caused by resid- 
uals of the compensation of the chromatic piston effect, which 
was most noticeable on the baseline A0-G1, and by systematic 
calibration uncertainties due to varying atmospheric conditions. 
Within the obtained errors of the UD diameter of about 4%, we 
do not see indications of any elliptical intensity distribution of 
f And. However, the ellipticity of ~4% expected for the night of 
September 18 is consistent with our data. 

Correction factors between UD diameter and limb-darkened 
(LD) disk diameters were computed using ATLAS 9 model at- 
mospheres (Kurucz 1993). For the spectral types of our target 
stars £ And and yU Peg and the wavelength range used for our ob- 
servations, we obtain values for ^ud/^ld of 0.974 and 0.976, re- 
spectively. The resulting LD diameters are #ld = 2.55+0.09 mas 
and 6» LD = 2.49 + 0.11 mas, respectively. Cohen et al. (119991 1 
give a diameter of 2.72 + 0.036 mas for £ And, based on spectro- 
photometry. This diameter is significantly larger, but the error 
smaller, than what was obtained in this work. Still, the spectro- 
photometric observations could be affected by the significant 
magnetic activity exhibited by ( And. The LD diameter of fi Peg 
obtained here is consistent with the earlier interferometric mea- 
surements of #ld = 2.53 ± 0.09 mas obtained with the NPOI and 
6>ld = 2.49+0.04 mas obtained with the Marklll interferometers 
(Nordgren et al. l2001l and Mozurkewich et al. 120031 1. increasing 
the confidence in the results presented here. 



4. Fundamental parameters 

4.1. Radius 



measured differential phase of the science target. We obtained 
results well within the adopted error. 

As a final data reduction step, the squared visibility ampli- 
tudes were calibrated for the interferometric transfer function, 
which was estimated using an average of the computed transfer 
functions based on the closest calibration star measurement be- 
fore that of each science target and the closest thereafter. The fi- 
nal error of the calibrated data includes the statistical error of the 
frames, the error in the correction for chromatic piston, and the 
standard deviation of the two transfer function measurements. 



3.2. Analysis of interferometric data 

Figure [T] shows the resulting visibility data of f And and of the 
check star Peg obtained from all three observing nights and 
compared to models of a uniform disk (UD). Because of the rel- 
ative large errors in the observations, no differences were seen 
in the measurements from different baselines. Thus all the base- 
lines were used together in the analysis. Table[3]lists the resulting 
uniform disk diameters of £ And and fi Peg for each of the nights 
separately, as well as for all three observing nights. When using 
the data from all the nights together, the diameter is estimated 
from the data obtained during the nights starting on September 
14 and 16, but the error is from all the data, i.e., including the 
data from the night starting on September 18. This is done be- 



The limb darkened diameter of £ And, obtained from the inter- 
ferometric observations, is 2.55 ± 0.09 mas. Together with the 
Hipparcos parallax of 17.24 ± 0.26 mas (van Leeuwen 2007) 
this can be used to determine the stellar radius with the follow- 
ing formula: R = ©ld^-, where Old is the limb-darkened angu- 
lar diameter in radians, 7r p the parallax in arcseconds, and C the 
conversion from parsecs to meters. For £ And this gives stellar 
radius of R = 15.9 + 0.8R Q , which is consistent with the 16.0 R Q 
estimated in Paper 1. 



4.2. Effective temperature 

The effective temperature of a star can be calculated from the 
interferometric diameter determination when combined with a 
bolometric flux measurement using the formula 



T e g - 



( 4/ b( 



1/4 



(1) 



where /bd is the bolometric flux and <x is the Stefan-Boltzmann 
constant. 

The bolometric flux of ( And was estimated using mea- 
surements on all the available photometric passbands and the 
getCal tool of the NASA Exoplanet Science Institute's interfer- 
ometric observation planning tool suite. The bolometric flux of 
/boi = 9.863 + 0.54 x 10~ 10 W/m 2 was obtained. Inserting this 
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Fig. 2. Doppler imaging results of £ And obtained with the TempMap^ code. In each plot the temperature map for an individual 
spectral line is shown with the observed spectra and photometry, including the corresponding fits to the data. In the middle panels, 
under the continuous line fits, the tiny vertical dashes are the l-cr error bars of the spectroscopic observations. 



value to Eq. Q] together with the limb-darkened angular diame- 
ter, gives T e s of 4665 + 140 K. This value is very close to, and 
the same within errors, as the r e ff « 4600 K used for Doppler 
imaging in Paper 1 and the current work. 



5. Doppler imaging 

For Doppler imaging we used the code TempMap, which was 
originally written by Rice et al. (1989). The code performs a 
full LTE spectrum synthesis by solving the equation of transfer 
through a set of ATLAS-9 (Kurucz [1993l) model atmospheres at 
all aspect angles and for a given set of chemical abundances. 
Simultaneous inversions of the spectral lines, as well as of 
the two photometric bandpasses, are then carried out using a 
maximum-entropy regularization. For the non-spherical ( And, 
a new version of the code was applied: TempMap £ (see Paper 1 
and the references therein) takes the distorted geometry of the 



evolved component in a close binary into account through the 
distortion parameter e. The elliptical distortion is approximated 
by a rotation ellipsoid, elongated towards the secondary star: 

e = -^1 - where a and b are the long and the short axes of 
the ellipsoid, respectively. The appropriate value of e for f And, 
0.27, as well as the overall system and stellar parameters, were 
adopted from Paper 1 (Table 2 therein). 

The 30 available UVES spectra (three exposures per night) 
covered 18 days, i.e., one full rotation cycle, thus allowing one 
Doppler reconstruction. The three nightly observations were av- 
eraged, since they were taken within approximately 120 sec. 
Thus, for further investigation we used the ten averaged spec- 
tra with an enhanced S/N value of ~600 or more (see Table|2]for 
more details on observations). 

Doppler imaging was performed for the well-known map- 
ping lines within the 6392-6440 A spectral range. Doppler maps 
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Fig. 3. Average map produced from all the maps shown in Fig. [2] 

for Fe i 6393, 6400, 641 1, 6421, 6430, and Ca i 6439 are shown 
in Fig |2] The individual maps revealed similar spot distribu- 
tions, i.e., mainly cool spots at low latitudes with temperature 
contrasts of 600-900 K with respect to the unspotted surface of 
4600 K. Cool polar features are also recovered, however, with 
significantly weaker contrast ranging from ~100K (Fei 6400) 
to a maximum of ~700K (cf. the Fei 6430 map). Numerous 
bright features also appear in the iron maps; however, as they 
occur near dominant cool spots they can be artifacts, so-called 
"rebound" features (see, e.g., Rice 2002]!. 

Despite the small difference between the temperature con- 
trasts of the respective maps and the spurious bright features, 
the resulting six Doppler maps are in very good agreement. This 
similarity is more conspicuous in Fig[3] where the average of the 
six individual maps is plotted. Averaging did not blur the overall 
structure. The most prominent feature is the belt of cool spots at 
the equatorial region, with the strongest concentration of spots 
located at the phase <p =0.75 and at another cool region ranging 
between phases <p = 0.0 - 0.4. Also a weak polar feature can 
be detected. This result is reminiscent of the result in Paper 1, 
where low-latitude dominant features also tended to concentrate 
at quadrature positions of opposite hemispheres for both observ- 
ing seasons. 

6. Discussion 

6. 1 . Comparison between spherical and ellipsoidal surface 
geometry in the inversion 

Another Doppler imaging code, INVERS7PD, which was writ- 
ten by N. Piskunov (see, e.g., Piskunov et al. Q990) and mod- 
ified by T. Hackman (Hackman et al. 2001), was also used to 
obtain a temperature map of ( And. In this inversion spherical 
geometry and only the Fe i 6400 A line were used. The obser- 
vations are compared to a grid of local line profiles calculated 
with the SPECTRUM spectral synthesis code (Gray & Corbally 
1994) and Kurucz model atmospheres (Kurucz 1993). In the cal- 
culations, 10 limb angles and nine temperatures between 3500 K 
and 5500 K were used. Photometry was not used as a constraint 
in this inversion as the ellipticity effect seen in the light curves 
cannot be properly taken into account when using spherical ge- 
ometry. 

The result of this inversion is shown in Fig. [4] The main spot 
structures and the temperature range are similar to the ones in 
the map obtained with TempMaPj: from the Fe 1 6400 A line. The 
cool spots concentrate on the equatorial region and especially on 
the quadrature points. The main spot is seen at the phase 0=0.75 
in the equatorial region, and there is also a prominent spot at 
the phase 0=0. 25 at higher latitudes. This is missing from the 
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Fig. 4. Temperature map of £ And obtained with INVERS7PD 
inversion code using spherical geometry and Fe i 6400 A line. 
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Fig. 5. Observed V magnitude of ( And (crosses) compared to 
the one calculated from the temperature map obtained using 
INVERS7PD (line). 

map obtained using TempMap £ , so it is most likely an artifact 
caused by using spherical geometry on an ellipsoidal star (cf. 
Fig. 3 in Paper 1). Furthermore, at the phase 0=0.25, the equa- 
torial region has a temperature close to that of the unspotted 
surface, unlike in the map obtained using ellipsoidal geometry. 
Also, the whole temperature scale is shifted by 100 K towards 
the cooler temperatures. On the whole, the temperature map ob- 
tained with the spherical geometry is very similar to the one ob- 
tained using ellipsoidal geometry and TempMap,.. As expected, 
the main differences occur at the quadrature points and espe- 
cially at the phase 0=0.25. Also, one has to keep in mind that 
the tests with TempMap e show that neglecting the ellipticity in 
the Doppler imaging reconstruction yields ~50-240% higher^ 2 
values in comparison to using the correct surface geometry. 

The main differences between the results from the spherical 
and ellipsoidal codes can be seen in the photometry. Figure [5] 
shows the normalised V observations of £ And compared to 
the ones calculated from the temperature map obtained with the 
code using spherical geometry. Only around the phases 0.4-0.6 
the two light curves show similar behaviour, and the photom- 
etry calculated from the INVERS7PD temperature map shows 
completely different behaviour than the observed one especially 
around the quadrature points. 

6.2. Chromospheric activity 

Chromospheric activity of £ And was investigated using the Ha 
line profiles, which appeared in absorption during the observa- 
tions, similarly to the other Balmer lines. Variations in the Ha 
line through the rotation cycle are shown in Fig. [6^. Both the red 
and the blue wings show strong variation at one, but different, 
phase. Also, most line profiles clearly show variable behaviour 
between velocities -100 km/s and +100 km/s. These variations 
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Fig. 6. Variations in the Ha profiles of £ And. a) The contin- 
uum normalised Ho- profiles, b) Residual profiles after the mean 
profile has been subtracted. The x-axis gives the velocity in com- 
parison to the mass centre of the giant component. 



are clearly seen already in the spectra, which have not been cor- 
rected to the continuum level. All the spectra show identical con- 
tinuum shapes, except approximately +5 A from the Ha line, cor- 
responding to the variation also seen in the normalised spectra 
used in the following analysis. 

To investigate the line-profile variations in more detail, the 
average profile was subtracted from all the profiles, thus creat- 
ing the residual profiles shown in Fig. [6j). The temporal vari- 
ations are clearly seen in these profiles. The most prominent 
features are the two strong absorption features seen at the ve- 
locities -350 - km/s and - 100 km/s. A dynamic spectrum, 
shown in Fig. [7] was also created from the difference profiles. 
Brighter colours in the plot correspond to enhanced emission 
and the darker colours to the enhanced absorption. The helio- 
centric Julian dates of the observations are shown with crosses 
in the plot. The data for the times where there are no observations 
are interpolations between the closest timepoints with data. The 
plotting over the heliocentric Julian date instead of the rotational 
phase was chosen, as some events are short lived, and the obser- 
vations in any case cover only slightly more than one rotation. 
The observational phases are given on the left side of the plot. 

The most noticeable feature in the dynamical spectrum is the 
strong enhanced absorption around the phase 1 .6 at the velocity 
+ 100 km/s. If the velocity seen in this absorption system was 
caused by the stellar rotation, it would be outside the stellar disk 
and would not be seen in absorption. Thus it must be a cloud 
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Fig. 7. Dynamic spectrum of the Ha line based on the residual 
profiles after the subtraction of the average line profile from the 
individual profiles. The x-axis gives the velocity and the y-axis 
the heliocentric Julian date on the left and phase on the right side. 
The crosses on the lefthand side of the plot give the heliocentric 
Julian dates of the observations and the dotted line the zero ve- 
locity. The brighter the colour, the more emission is observed. 



of cool gas in the stellar atmosphere that is falling into the star. 
This could be the final stages of a flare event. No evidence of 
such an event is seen in the earlier observations, but it could have 
occurred during the one-night gap in the observations. More en- 
hanced absorption is seen at the first observation (phase 0.94) ex- 
tending to the very blue, to -300 km/s and beyond. This could be 
caused by a mass ejection event with a strong line of sight com- 
ponent. In the following observation, more enhanced absorption 
is seen spanning the velocities -100 - +100 km/s. 

Enhanced emission occurs at three main locations: around 
the phases 1 .15-1 .30 at the velocities -40 - -70 km/s, around the 
phases 1 .2-1 .4 at the velocities +60 - 100 km/s and at the phases 
1.7-1.8 at the velocities +40 - 70 km/s. These features have ve- 
locities that place them slightly outside the stellar disk, and thus 
they could be caused by prominences seen at the stellar limb. 
The prominence seen at the blue edge around phases 1.1-1.3 
is most likely the same one as seen at the red edge 0.5 in a later 
phase (i.e., at phases 1.7-1.8). Also, a weak enhanced absorption 
feature is seen at phase 1.4 around the velocity -20 km/s, which 
could be caused by the prominence starting to cross the stellar 
disk. This prominence could be centred around phase 1 .5, which 
in the binary reference frame is the phase pointing away from the 
secondary. The enhanced emission seen in the red around phases 
1.2-1.4 are, based on their velocities, also most likely caused by 
prominences. However, they have to be short lived in nature, as 
no evidence of them is seen in the observations before or after. 
These prominences would be at the disk centre approximately at 
phases 1.0 and 1.1, which places them on the side phasing the 
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Fig. 8. All the differential V and y magnitudes of £ And in the 
instrumental system plotted against the rotational phase. 
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Fig. 9. All the differential V and y magnitudes plotted against 
the heliocentric Julian date. The large crosses show the mean 
magnitudes for that time. 



secondary. They also coincide with the weaker cool region seen 
around the phases 0.0-0.4 in the Doppler image. 

6.3. Long-term magnetic activity of % And 

The long-term activity in ( And is investigated based on the pho- 
tometric V and y band observations obtained with the Wolfgang 
and Amadeus automatic photometric telescopes. The observa- 
tions between December 1996 and October 2002 were already 
used in Paper 1. Here, we also use observations obtained be- 
tween June 27, 2003 and October 25, 2008, in total 211 new V 
magnitudes. 

When all the instrumental differential magnitudes are plot- 
ted against the phase, see plot Fig. [8] the variation caused by 
the ellipticity effect is clearly seen. Still, the observations show 
much larger scatter around the ellipticity curve than is expected 
from the measurement error of 0.01-0.02 magnitudes. This in- 
dicates that there are also significant variations due to starspots. 
Evidence of changes in the activity level are also seen when all 
the observations are plotted against the Julian Date in Fig. [9] In 
this plot the small crosses give the individual observations and 
the large crosses the mean of that time period. No mean is given 
for some time periods, as there are so few measurements, or they 
are grouped such, that the full light-curve was not sampled. 

Changes in the mean magnitudes could be interpreted as a 
solar-like activity cycle. A spectral analysis of the mean mea- 
surements was carried out using the Lomb method (Press et 
al. 1 1992| l. The result indicates a presence of a cycle with a cy- 
cle length of 5.9 + 0.7 years, but the false alarm probability is 



0.36. Thus the period can be spurious, and more measurements 
are needed for confirming it. 

As can be seen from Fig. [9] £ And is on average brighter by 
0.014 magnitudes during the VLT observations presented here 
than during the December 1997 - January 1998 KPNO obser- 
vations presented in Paper 1. This implies that the spot cover- 
age, and/or the spot temperature, is different between the two 
epochs. For studying this further, the temperature maps obtained 
with the Ca i 6439 A line, were investigated both from the VLT 
and KPNO observations. The hottest temperature is basically the 
same in both maps, 4600 K for VLT and 4620 K for KPNO. Still, 
the coolest temperatures are very different. For the VLT obser- 
vations, the coolest spots are 3710 K and for the KNPO map 
3480 K. Also, the number of surface elements having a tem- 
perature of 4000 K or less in the map obtained from the VLT 
is 75% of the surface elements with those temperatures in the 
KPNO map. This investigation supports the existence of a activ- 
ity cycle, which is also indicated by the long-term photometry 
of £ And. Still, one must keep in mind that the temperatures in 
the Doppler images are very sensitive to the data quality, and the 
VLT data are superior to the KPNO ones. 



7. Conclusions 

The following conclusions can be drawn from the optical in- 
terferometry, high-resolution spectroscopy and broad band pho- 
tometry presented in this work. 

1. Optical interferometry gives an apparent diameter of 2.55 + 
0.09 mas for £ And. Using the Hipparcos parallax, this trans- 
lates into a stellar radius of R = 15.9 ± 0.8R Q , which is in 
line with the earlier radius determinations. 

2. Combining the interferometrically determined diameter and 
bolometric flux gives an effective temperature of r e jf = 
4665 + 140, which is consistent with the values determined 
through Doppler imaging. 

3. The expected ellipsoidal stellar geometry with ~4% differ- 
ence between the long and short axes cannot be confirmed 
with the current interferometric observations, which have er- 
rors of about 4% in the diameter measurement. However, the 
highest ellipticity expected for the night of September 18 is 
consistent with the data. 

4. The Doppler images reveal cool spots on the surface of the 
primary of the £ And binary. The spots are located in the 
equatorial region, and the main concentration of spots is seen 
around phase 0.75, i.e., 0.25 in phase from the secondary. 
Another weaker cool region spans the phases 0.0-0.4, again 
around the equator. There are also indications of a cool polar 
cap. On the whole, this spot configuration is very similar to 
the one seen in the earlier published 1997/1998 data. 

5. Long-term photometric observations indicate an activity cy- 
cle, but more measurements are needed to confirm this and 
its period. The investigation of the Doppler maps obtained 
January 1998 and September 2008 also hint at an activity 
cycle. 

6. The chromospheric activity, investigated from the Ha-line, 
shows evidence of both prominences and cool clouds. The 
prominences do not seem to show any strong evidence of oc- 
curring at certain locations in the binary reference frame, nor 
are they associated with the coolest spot seen on the surface. 
On the other hand, one of the detected prominences seems to 
be related to the group of weaker cool spots located at phases 
0.0-0.4. 
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Table 1. Log of the VLTI/ AMBER observations. 



Date 


Time 


Target 


Purpose 


FINITO 


AM 


Seeing 


TO 


Sep. 


UTC 










["] 


[msec] 


15 


4:24-4:30 


<?Psc 


Calibrator 


on 


1.17 


0.66 


3.1 


15 


4:47-4:51 


SPsc 


Calibrator 


off 


1.17 


0.66 


3.0 


15 


5:04-5:07 


/"Peg 


Check star 


on 


1.63 


0.65 


3.1 


15 


5:12-5:16 


/"Peg 


Check star 


off 


1.65 


0.62 


3.2 


15 


5:24-5:27 


41Psc 


Calibrator 


on 


1.19 


0.67 


2.9 


15 


5:32-5:35 


41Psc 


Calibrator 


off 


1.19 


0.69 


2.8 


15 


5:45-5:49 


/"Peg 


Check star 


on 


1.79 


0.62 


3.1 


15 


5:55-5:58 


/"Peg 


Check star 


off 


1.85 


0.54 


3.6 


15 


6:03-6:06 


/"Peg 


Check star 


on 


1.90 


0.69 


2.8 


15 


6:14-6:17 


(9Psc 


Calibrator 


on 


1.30 


0.71 


2.7 


15 


6:22-6:25 


0Psc 


Calibrator 


off 


1.33 


0.62 


3.1 


15 


6:46-6:49 


( And 


Science target 


on 


1.58 


0.71 


2.9 


15 


6:53-6:57 


£And 


Science target 


off 


1.60 


0.79 


2.6 


15 


7:02-7:06 


41Psc 


Calibrator 


on 


1.32 


0.71 


2.9 


15 


7:13-7:22 


£And 


Science target 


on 


1.66 


0.65 


3.2 


15 


7:29-7:33 


HD7087 


Calibrator 


on 


1.53 


0.80 


3.2 


15 


7:38-7:41 


HD7087 


Calibrator 


off 


1.55 


0.82 


3.2 


15 


7:53-7:56 


£And 


Science target 


on 


1.84 


0.97 


2.2 


15 


8:05-8:08 


HD 15694 


Calibrator 


on 


1.13 


0.73 


2.9 


17 


5:11-5:14 


0Psc 


Calibrator 


on 


1.19 


0.91 


2.4 


17 


5:42-5:46 


/"Peg 


Check star 


on 


1.82 


1.02 


2.1 


17 


6:08-6:15 


41Psc 


Calibrator 


on 


1.22 


1.09 


2.0 


17 


6:36-6:40 


/"Peg 


Check star 


on 


2.30 


0.97 


2.5 


17 


6:56-6:57 


4" And 


Science target 


on 


1.62 


0.85 


2.6 


17 


7:09-7:13 


HD7087 


Calibrator 


on 


1.50 


0.86 


2.5 


17 


7:25-7:27 


£And 


Science target 


on 


1.74 


0.86 


2.5 


17 


7:41-7:44 


HD 15694 


Calibrator 


on 


1.12 


0.91 


2.4 


17 


7:57-8:00 


£And 


Science target 


on 


1.92 


0.73 


2.9 


19 


3:47-3:53 


(9Psc 


Calibrator 


on 


1.18 


0.74 


1.9 


19 


5:13-5:16 


/"Peg 


Check star 


on 


1.72 


0.78 


1.8 


19 


5:28-5:31 


0Psc 


Calibrator 


on 


1.23 


0.75 


1.8 


19 


5:46-5:49 


/"Peg 


Check star 


on 


1.90 


1.06 


1.3 


19 


6:01-6:04 


41Psc 


Calibrator 


on 


1.22 


0.97 


1.4 


19 


6:19-6:21 


/"Peg 


Check star 


on 


2.18 


0.89 


1.3 


19 


6:54-6:57 


<TAnd 


Science target 


on 


1.64 


0.81 


1.7 


19 


7:19-7:22 


HD 7087 


Calibrator 


on 


1.54 


1.18 


1.2 


19 


7:52-7:55 


(And 


Science target 


on 


1.95 


1.55 


0.9 


19 


8:18-8:22 


HD 15694 


Calibrator 


on 


1.16 


1.64 


0.9 



Table 2. The high-resolution spectroscopy with UVES at VLT. 



Date 


HJD 


phase 


S/N 




2450000+ 






13.09.2008 


4722.81672 


0.940 


866 


15.09.2008 


4724.70463 


0.046 


650 


17.09.2008 


4726.71012 


0.159 


586 


18.09.2008 


4727.68700 


0.214 


628 


19.09.2008 


4728.73866 


0.273 


762 


21.09.2008 


4730.72025 


0.384 


693 


23.09.2008 


4732.72873 


0.497 


914 


25.09.2008 


4734.65166 


0.606 


785 


27.09.2008 


4736.69360 


0.721 


612 


01.10.2008 


4740.75333 


0.949 


802 



